Effect of the Coulomb repulsion on the ac transport through a quantum dot 
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We calculate in a linear response the admittance of a quantum dot out of equilibrium. The 
interaction between two electrons with opposite spins simultaneously residing on the resonant level 
is modeled by an Anderson Hamiltonian. The electron correlations lead to the appearence of a new 
feature in the frequency dependence of the conductance. For certain parameter values there are two 
crossover frequencies between a capacitive and an inductive behavior of the imaginary part of the 
admittance. The experimental implications of the obtained results are briefly discussed. 
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The double-barrier resonant tunneling systems (DRBTS) are in the focus of intensive experimental [Q and theo- 
£>v retical [||-|| investigations. While the main characteristics of the steady-state transport properties of such systems 
have been well inderstood in terms of several basic approaches — the kinetic equation approach , by a Landauer- 
Biittiker-type formula ||, and via the Wigner- function Q - some of the features in the time-dependent behavior of 
the DRBTS still remain unclear. In particular, the long lasting question of the fact whether the tunneling transport 
, through the structure is via an establishing a coherent wave function, or can be viewed as a sequence of quantum 
tunneling events flC| , |ll|| . The answer to this question is of major importance to the practical applicability of these 
devices as a high frequency resonator s 02 -15 
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■^j- ■ In a recent paper Fu and Dudley |14|| Tiave utilized a model of non-interacting electrons transmitted through the 
t-H ' structure under the influence of a small ac bias superimposed upon the driving dc voltage, having shown that for 
some values of the system's parameters the frequency-driven behavior can be simulated by an equivalent electrical 
circuit with an additional inductive element - feature, which was earlier encountered in the modelling the frequency 
dependence of the DRBTS ( Jl(],[ll]] ) . In this context, however, it is important to correctly account for the Coulomb 
interaction effects due to the repulsion between two electrons simultaneously residing on the resonant level - phe- 
nomenon, which was shown to give rise of a numerous effects in the dc transport ^U6|J^]. In a previous paper [ fl5| we 
have calculated the linear response admittance of a quantum dot with interacting electrons in the case of a zero dc 
voltage, thus isolating the effects due to the Coulomb on-site repusion in an equilibrium state of the system. We have 
obtained a feature in both the real and the imaginary part of the admittance at £1 ~ E c . In a recent paper Brandes, 
Weinmann, and Kramer |l6| have studied the ac conductance of a tunnel junction in a linear response approximation. 
Having introduced the picture of a simultaneous transmission of a electron-hole pair through the system in the ac 
transport regime, they argued that at high frequencies the Coulomb interaction would be of minor importance since 
the pair is electro-neutral. They estimated the upper frequency limit for the Coulomb interaction to be noticeable 
luq ks 1GHz for the currently accsessible experimental parameters. 

In this work we apply a non-equilibrium technique to this problem to deal with the case of a non-zero applied dc 
voltage. We study the effects of the Coulomb interaction on the relatively low-frequency behavior of the conductance 
through the dot. It is worth mentioning at this point that the formalism developed here is equally applicable either to 
the case of symmetrical (Tl = Tr) or asymmetrical (Tj, ^ Tr) coupling of the well to the leads [[17|, and can be easily 
extended to take into account additional effects due to environmental fluctuations (i.e. electron-phonon coupling). 
The conductance is calculated as a function of the dc voltage and the frequency f2 of the applied ac bias for f2 of the 
order of several 7, where 7 is the elastic width of the resonant level. 

We found that the Coulomb repulsion has a profound effect on both the frequency-driven conductance and the 
energy losses over a relatively broad range of parameters. In particular, a new feature in the ac response is found 
in the case of finite Coulomb repulsion energy E c , which we contribute to the fundamental way in which the well- 
electrons energy spectrum cnanges in the presence of electron-electron interactions - and this is in marked contrast 
to the non-interacting (E c — > 0) limit. The imaginary part of the admittance behaves in a way consistent with the 
Kramers-Kronig dispersion relations, thus giving rise to an additional featute for a finite E c . We further discuss in 
some details the experimental implications of the theoretical findings reported here, giving the conditions necessary 
for observing the predicted effects, and discussing some recent experimental results as well. 
The Hamiltonian of the quantum well, coupled to the leads, can be written as 

H = H e £ a L afc - + e P b U b p° + e -Yl c ° c ° + E ^ n i + Y,( TLkC ° ak - + h - c ) + Y,^ TR p b U c - + h - c ~)- w 

kcr per a kcr per 

It is expressed in terms of the creation (annihilation) operators in the emitter a^ a , collector b p(7 , and the well (c CT ) with 
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k(p) the corresponding quasimomenta and a =t , | the spin index. Here and are the single-particle energies in 
the emitter and the collector leads, respectively. E c is the Coulomb repulsion between electrons with opposite spins. 
The single-electron energies are measured from the corresponding Fermi levels and in the emitter and collector, 
and the dc bias is /i£ — = eV. e c = ej? + aeV, where t° c is the bare resonant level energy and a measures the portion 
of the voltage drop on the quantum well(a ~ 0.5 for the symmetric structure T^k — Tr p ). The particle- number 
operator is n a — c\c a . T^f. and Tr p are the tunneling matrix elements which are bias dependent. 

In order to investigate the response of the system to an external perturbation we introduce an additional ac voltage, 
which is superimposed on the dc bias eV. The Hamiltonian of this ac signal has the form: 

H ex t(t) = -aeu(t) ^ 4 c <r - eu(t) ^ b \a b P<y + h - c -, ( 2 ) 

a pa 

where the field operators are taken at the time t. The external field is u(t) ~ u^e and its amplitude is assumed to 
be very small in order not to disturb the electron distribution in the well. In writing this form of the H ext we have 
set fi R = (13). 

First we calculate the electron Green's functions. We apply the non-equilibrium (Keldysh) formalism where one 
introduces the retarded(advanced) and distribution Green's functions [|l8|,[l9| . We assume that the relaxation processes 
in the leads are much faster than in the quantum well. Thus, we can consider the leads as equilibrium systems and 
the corresponding Green's functions are given by the usual expressions for a non-interacting electron system 1 18 lj| . 

The retarded (advanced) G r r a ^(uj) and the distribution G < (ui) Green's function for the electrons in the well corre- 
sponding to the full Hamiltonian H are determined using the irreducible Green's function method [H . This method 
has been applied successfully to the Hubbard model |2lj . Its main advantage consists of treating all the truncations 
of the higher order Green's functions with the same accuracy - consistent with already choosen algebra of relevant to 
the problem operators. It is worth mentioning that it is valid both in the weak and the strong correlation limits. 

We recall that in the non-equilibrium formalism the Green's functions depend on two time variables ti,t%. It can 
be shown that the retarded (advanced) Green's functions depend only on the "relative" time t = t\ —ti while the 
distribution Green's function depends on both t and T — ll ^ t2 |^9[ . 

The well electrons Green's function is obtained in the form (for more details see pj|) 

G r (uj) = — — \ — -. (3) 

(cj - e c - Eo)(w - e c - So - Si) - E c (u> - e c - S - £2) 

The explicit expressions for the self-energy parts So, Si, S2 in Eqn. (|J) can be found in Ref. 15. We use the following 
notation: e ca = e c + E c (l— < n_ CT >) where < n a > is the average number of well electrons with spin a. 

This Green's function describes two energy levels for the quantum well electrons - a lower level with energy e c and 
an upper level with energy E c + e c . 

We should stress that the derivation of the Green's function presented here is valid for temperatures higher than the 
characteristic temperature for this problem - the Kondo temperature Tk- Lacroix p2| has shown that for temperatures 
T < Tk this truncation procedure omits terms which are divergent at the Fermi level. These terms give rise to the 
Kondo effect. 

The distribution Green's function is calculated assuming that all transient processes after the switching on the dc 
bias have decayed. In this case one can consider the distribution Green's function to be independent on the time T. 
The result for G < (lu) is found in the following closed form: 

G < (u)=-F(u){G r (v)-G a (w)) (4) 
where 

J2\T Lk \ 2 A^k^) + J2\Tr p \ 2 B^u) 
~ E \T Lk \ 2 (A r (k, lj) - A a (k, lj)) + E |T flp | 2 (B r (p, to) - B a (p, u )) 

k p 

is the new non-equilibrium (but steady-state) distribution of electrons in the quantum well. In this expression 
A r ( a ) (B r ( a )) an d A< (£?<) are the retarded (advanced) and the distribution Green's functions in the left (right) lead, 
respectively. 

Next we calculate the current through the quantum well in the presence of a time-dependent electric field. The 
total current is given by I = (II + Ir)/2 where Il(r) is the current through the left(right) barrier, respectively 
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I(t) = [T Lk cUt)a k „(t)-T£ k al(t) C(r (t) 

k(7 

+ E [TR P bl(t)cAt) - T* Rp cUt)b pa (t)]). (6) 

pa 

To obtain this expression we have used the Shockley - Ramo theorem p3f| . 

In order to calculate the current i(t) due to the external ac voltage u(t), we use the linear response formalism. In 
this way we obtain the admittance of the quantum dot Y(u>) = i(u>)/u(u>) where u(lj) = uq(8(lo — 57) + 5(lu + SI)) is 
the Fourier transform of the external ac bias. The explicit expression for the admittance is obtained in the form: 

F(O) = J dto {~2i (a + D r (Lo, SI)) ( lL (cu)f L (io) - *y R (u>)f R (w)) G r (u + (l)G r (u) 

+ 2i(a + D a (u, SI)) ( 7L (w + Q)f L (u + SI) - j r (lj + Q)f R (u + SI)) G a (uj + Q)G a (u) 

- [{a + D a {uj + SI)) F(uj + SI) (G r (w + Si) - G a {uj + SI)) G a (u) (7) 
+ (a + D r (ui + SI)) F{uj)G r {uj + Si) (G r (cu) - G a (u))] 

x (A r (u + St) - A a (uj) - B r {u + Sl) + B a (u>))} (8) 
In writing Eqn. Q we have used the following notations: 
A r{a) (cj) = ^|T Lfc | 2 A r(a) (fc, w ), (9) 



B r (a){u) = ^2 \T Rp \ 2 B r{a) (p,uj), (10) 
p 

D r ( a )(v,Sl) =^2\T Rp \ 2 B r{a) (p,ui + Sl)B r(a) (p,uj). (11) 
p 

The tunneling matrices and T Rp are related to the level widths for the leads' electrons via the usual equations 
7l(w) = 7r \TLk\ 2 fii^ — ^k), 7fl(w) = 7T \T Rp \ 2 5(u! — e p ). In the following 71,(0;), ~f R (u)) are taken to be independent 

k p 

of u>. 

Now we present our numerical results for the ac conductance a(Sl) = ReY(Sl) and the energy losses ImY(Si) through 
the quantum well. We calculate them for T > Tk assuming that E c ^> T, 7 = 7^ +j R and T 3> 7. First we solve self- 
consistently the equation for the average number of quantum well electrons (n a ) = {n~ a ) = n = — J du>/2irImG < (uj) 
(we consider a non-magnetic solution). We take a broad flat density of states for the leads' electrons. In Fig. 1 
we show the dynamical conductance cr(0, V) calculated for a bare level energy e' '' = 0.2E C and for a symmetrical 
coupling of the leads to the well - 7l = 7h- 

For O — > (dc limit) cr(V) has two maximums. This structure reflects the energy spectrum of the well electrons 
- there are two channels for the electrons to tunnel through the well. When the number of well electrons is smaller 
than one (n < 0.5) the tunneling is predominantly through the lower (resonant) level. For n > 0.5 the lower level is 
filled with electrons and they are transfered through the upper level. 

In this paper we show that the effect of the electron correlations (the Coulomb repulsion E c ) can also be observed 
in the SI dependence of the conductance. For relativelylow frequency (SI ~ 2 — 37) the conductance decreases with 
Si (similarly to the case of non-interacting electrons JlJ,|lJ]) since the electrons cannot follow the applied ac field. 
The electrons build-up in the well and fill the upper level (there is more than one electron in the well). This opens 
an additional tunneling channel through the upper level. Consequently for higher frequencies the tunneling current 
increases (the feature at ~ 8 — IO7 for relatively low dc voltage). For larger frequencies the conductance again falls 
off with SI, it becomes negative and tends to zero with negative values. This behaviour is in marked contrast to the 
non-interacting case. In the latter the conductance is positive (for a resonant level above the right chemical potential) 
and never changes sign. 

When the dc voltage increases the feature in the SI dependence of a we have just discussed is almost smeared out. 
The conductance is a monotonically decreasing function of SI but the non-zero response to the ac field spreads to 
substantially higher frequencies compared to the non-interacting case. 
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For dc voltages in the region of the second peak (at SI — > 0) when the renormalized level is brought well above hr 
the conductance is a non-monotonic function of SI. It slightly increases for SI < 57 and then decreases with non-zero 
values of <r(Sl) up to SI ~ 157. 

In Fig. 2 we show the imaginary part of the admittance Im Y(£l) calculated for two typical cases: a) — —0.1 12E C 
and b) = — 0.01E C . The applied dc voltage was taken to be eV = 0.2E C and 7 = 0.04E C . When the upper level 
is above the right chemical potential the admittance shows a capacitive behavior - Im Y(Q) > 0( Fig. 2, curve a)). 
Note that the sign of the imaginary part of the admittance is opposite to the sign in Ref . p4[ . 

Particularly interesting is the result presented in Fig. 2, curve b). In this case \e c — (j,r\ < 7/2. It shows that there 
are two frequencies at which a crossover from a capacitive to an inductive behavior (and vice versa) is obtained. At 
low frequencies the admittance is inductive, then it changes sign and this is the crossover to a capacitive admittance. 
Fu and Dudley [ fl4[ studied the same case in the non-interacting picture of the resonant tunneling. They obtained an 
inductive behavior - Im Y(fl) < and no crossover to a capacitive one. Moreover, in the interacting case there is one 
more crossover frequency at which the admittance changes back from a capacitive to an inductive. When V = this 
frequency is SI ~ ec + E c — (ir. With increasing V the first crossover frequency diminishes and the second increases 
so that at sufficiently high voltage the behaviour of the admittance will be inductive. The appearence of the second 
crossover frequency is in agreement with the conclusion of Brandes, Weinmann, and Kramer Jl6[ ] - at high frequencies 
the behavior of the admittance is not strongly affected by the electron interactions except the feature at SI ~ E c 
discussed in [ p^[ . 

When |e c — /ir\ > 7/2 (as in Fig. 2, curve a) Fu and Dudley obtained a crossover to an inductive admittance 
at SI ~ |e c — With E c ^ this is possible only if e c is well below ^r and the upper level is above ^jlr 

For a quantum dot with both the resonant and the upper level well below the right chemical potential for low V the 
admittance is inductive. A crossover to a capacitive behavior can be obtained for sufficiently high V at SI ~ e c +E c — [lr 
(then the upper level will be moved above fiR) |24}] . 

All these results clearly indicate that as in the non-interacting case JlJ] the frequency behavior of DBRTS cannot 
generally be simulated by any LRC equivalent-circuit model. 

To address the experimental detection of the the effects under consideration here, let us recall the range of the 
three main parameters involved in the model -£",,> T > 7. The on-site repulsion energy can be estimated from 
E c ~ e 2 /eL, where L is the size of the confined region, and e is the dielectric constant for the GaAs. Thus for a 
quantum dot with average size of 100 A, one gets E c ~ 1 — lOmeV^. It is less evident how to estimate the elastic 
broadening constant 7, but we can use the estimation given in Ref. pf| for 7 ~ 10 — 20[ieV for a structure of 
about the same size. Therefore, the type of effects discussed in this study to be detected experimentally, one needs 
temperatures of few degrees K and a external frequency range up to hundreds GHz - requirements accessible at the 
present time. Let us only mention at this place that because of the temperature ranges used (T = 77K, and room 
temperature, accordingly) the results of two recent ac experiments (see Refs. |l^,|ll|) did not show the frequency 
dependence discussed here. This is because the measurements were performed in the regime where the thermal 
fluctuations prevail over the "Coulomb blockade" effect ]26| ]. 

We should point out that the expression for the current Eqn. (JsJ) does not include the displacement currents through 
the parasitic capacitances. They can be accounted for by considering some electrostatic model of the dot. Our aim 
was to extract the effects of the electron correlations comparing our results with the available works on non-interacting 
electron tunneling where the displacement currents were not included. One can expect that the parasitic currents 
would modify the relatively high frequency behavior of the admittance |2^| . 

In conclusion, we have calculated the linear-response admittance of a quantum dot with interacting electrons. We 
show that the effect of the electron correlations can be observed in the frequency dependence of the conductance 
where a new feature appears when the electrons tunnel through the upper level. For a certain parameter values there 
are two frequencies at which the imaginary part of the admittance changes sign i.e. a crossover between a capacitive 
and an inductive behavior. We discuss the experimental conditions for observation of these effects. 
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FIG. 1. Conductance a as a function of the applied dc bias across the structure and the external frequency f2 calculated 
for a quantum dot with = 0.2E C . 



FIG. 2. The imaginary part of the dynamical admittance ImY as a function of the frequency Q of the external ac signal 
for two positions of the bare energy level: a) ei ' = — 0.112_E C (the solid line), and b) ei°^ = — 0.01-E C (the dashed line). The 
applied dc voltage is eV = 0.2E C , and 7 = 0.04i5 c . 
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